Kondo quantum dot coupled to ferromagnetic leads: 
a study by numerical renormalization group technique 
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We systematically study the influence of ferromagnetic leads on the Kondo resonance in a quan- 
tum dot tuned to the local moment regime. We employ Wilson's numerical renormalization group 
method, extended to handle leads with a spin asymmetric density of states, to identify the effects 
of (i) a finite spin polarization in the leads (at the Fermi-surface), (ii) a Stoner splitting in the 
bands (governed by the band edges) and (iii) an arbitrary shape of the leads density of states. For 
a generic lead density of states the quantum dot favors being occupied by a particular spin-species 
due to exchange interaction with ferromagnetic leads leading to a suppression and splitting of the 
Kondo resonance. The application of a magnetic field can compensate this asymmetry restoring 
the Kondo effect. We study both the gate-voltage dependence (for a fixed band structure in the 
leads) and the spin polarization dependence (for fixed gate voltage) of this compensation field for 
various types of bands. Interestingly, we find that the full recovery of the Kondo resonance of a 
quantum dot in presence of leads with an energy dependent density of states is not only possible by 
an appropriately tuned external magnetic field but also via an appropriately tuned gate voltage. For 
flat bands simple formulas for the splitting of the local level as a function of the spin polarization 
and gate voltage are given. 

PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 



I. INTRODUCTION 

The interplay between different many-body phenom- 
ena, such as superconductivity, ferromagnetism, or the 
Kondo effect, has recently attracted a lot of experimen- 
tal and theoretical attention. A recent experiment of 
Buitelaar 1 et al. nicely demonstrated that Kondo corre- 
lations compete with superconductivity in the leads. The 
interplay between Kondo correlations and itinerant elec- 
tron ferromagnetism in the electrodes, has theoretically 
been intensively studied within the last years, 2-9 initially 
leading to controversial conclusions. 

For effectively single-level quantum dots (i.e. dots with 
a level spacing much bigger than the level broadening T), 
consensus was found that a finite spin asymmetry in the 
density of states in the leads results (in general) in a 
splitting and suppression of the Kondo resonance. This 
is due the spin dependent broadening and renormaliza- 
tion of the dot level position induced by spin-dependent 
quantum charge fluctuations. In terms of the Kondo spin 
model it can be treated as an effective exchange interac- 
tion between a localized spin on the dot and ferromag- 
netic leads. Moreover, it was shown that a strong cou- 
pling Kondo fixed point with a reduced Kondo temper- 
ature can develop 6, even though the dot is coupled to 
ferromagnetic leads, given an external magnetic field 6 or 
electric field 7 (gate-voltage) is tuned appropriately. Ob- 



viously, in the limit of fully spin polarized leads, when 
only one spin component is present for energies close to 
Fermi surface (half- metallic leads), the effective screen- 
ing of the impurity cannot take place any more and the 
Kondo resonance does not develop. A part of these the- 
oretical predictions have recently been confirmed in an 
experiment by Pasupathy 10 et al. The presence of ferro- 
magnetic leads could also nicely explain the experimental 
findings of Nygard 11 et al. 

Since the interplay between ferromagnetism and strong 
correlation effects is one of the important issues in spin- 
tronics applications, there are currently many research 
activities going on in this direction. The goal is to ma- 
nipulate the magnetization of a local quantum dot (i.e. 
its local spin) by means of an external parameter, such as 
an external magnetic field or an electric field (a gate volt- 
age) , with a high accuracy. This would provide a possible 
method of writing information in a magnetic memory. 12 
Since it is extremely difficult to confine a magnetic field 
such that it only affects the quantum dot under study, 
it is of big importance to search for alternative possibil- 
ities for such a manipulation (e. g. by means of a local 
gate- voltage as proposed by the authors 13 ). 

In this paper we push forward our previous work by 
performing a systematic analysis on the dependence of 
physical quantities on different band structure proper- 
ties. Starting from the simplest case we add the ingre- 
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dients of a realistic model one-by-one allowing a deeper 
understanding of the interplay of Kondo model and itin- 
erant electron ferromagnetism. In this paper we ex- 
tend our recent studies carried out in this direction 6,7 ' 13 
and illustrate the strength of the analytical methods by 
comparing the results predicted by them to the results 
obtained by the exact numerical renormalization group 
(NRG) method. 14 While in Refs.[2-9] the dot was at- 
tached to ferromagnetic leads with an unrealistic, spin- 
independent and flat band - with a spin-dependent tun- 
neling amplitude - we generalize this treatment here by 
allowing for arbitrary density of states (DOS) shapes. 13 
In particular, we carefully analyze the consequences of 
typical DOS shapes in the leads on the Kondo resonance. 
We explain the difference between these shapes and pro- 
vide simple formulas (based on perturbative scaling anal- 
ysis 15,16 ) that explain the numerical results analytically. 

We study both the effect of a finite leads spin polar- 
ization and the gate voltage dependence of a single-level 
quantum dot contacted to ferromagnetic leads with three 
relevant DOS classes: (i) for flat bands without Stoner- 
splitting, (ii) for flat bands with Stoner-splitting and (iii) 
for an energy dependent DOS (also including Stoner- 
splitting). For this sake we employ an extended version 
of the NRG method to handle arbitrary shaped bands. 17 

The article is organized as follows: In Sec. II wc define 
the model Hamiltonian of the quantum dot coupled to 
ferromagnetic leads. In Sec. Ill we explain details of the 
Wilson's mapping on the semi-infinite chain in the case 
of the spin-dependent density of states with arbitrary en- 
ergy dependence. Using the perturbative scaling analysis 
we give prediction for a spin-splitting energy for various 
band shapes in Sec. IV. In Sec. V the results for spin- 
dependent flat DOS are demonstrated together with the 
Friedel sum rule analysis. The effect of the Stoner split- 
ting is discussed in Sec. VI together with comparison to 
experimental results from Ref.[ll] and an arbitrary band 
structure in Sec. VII. We summarize our findings then 
in Sec. VIII. 



II. MODEL: QUANTUM DOT COUPLED TO 
FERROMAGNETIC LEADS 

We model the problem at hand by means of a single- 
level dot of energy £d (tunable via an external gate- 
voltage Vq) and charging energy U that is coupled 
to identical, nonintcracting leads (in equilibrium) with 
Fermi-energy p — 0. Accordingly the system is described 
by the following Anderson impurity model 

H = He + Tied + Hd , 

H d = e d ^h a + Uh^hi- BS Z , (1) 

(7 

with the lead and the tunneling part of the Hamiltonian 

'Hi = y^^rkacikgCrka , (2) 
rkcr 



Hid = ^{Vrkdicrka + h.C.) . (3) 
rka 

Here c r ka and d a (n a = rfJ.rf CT ) are the Fermi operators 
for electrons with momentum k and spin a in lead r (r = 
L/R) and in the dot, respectively. The spin-dependent 
dispersion in lead r, parametrized by e r k<j, reflects the 
spin-dependent DOS, p r<J (uj) = ^ fe <5 (u> — e r ka), in lead 
r; all information about energy and spin dependency in 
lead r is contained in the dispersion function e r ka- V r k 
labels the tunneling matrix-element between the impu- 
rity and lead r, S z = (h^ — h{)/2, and the last term in 
Eq. (1) denotes the Zeeman energy due to external mag- 
netic field B acting on the dot spin only. Here we neglect 
the effect of an external magnetic field on the leads' elec- 
tronic structure as well as a stray magnetic field from 
the ferromagnetic leads. The coupling between the dot 
level and electrons in lead r leads to a broadening and 
a shift of the level £d, £d — » ?d (where the tilde denotes 
the renormalized level) . The energy and spin dependency 
of the broadening and the shift, determined by the cou- 
pling, T ra (uj) = np rrT (uj)\V r (ui)\ 2 , plays the key role in 
the effects outlined in this paper. Henceforth, we assume 
V r k to be real and k- independent, V r k a — V r , and lump 
all energy and spin-dependence of T r<7 (uj) into the DOS 
in lead r, p rcr (uj). 18 

Without loosing generality we assume the coupling 
to be symmetric, VL = Vr. Accordingly, by perform- 
ing a unitary transformation 19 He simplifies to Hi = 

\Hka ( e Lfe CT + eRk<r)al ka a skr r, where a sk a denotes the 
proper unitary combination of lead operators which cou- 
ple to the quantum dot and we dropped the part of the 
lead Hamiltonian which is decoupl ed from t he dot. With 
the help of the definitions V = \JVl + V^, aka = a s ka 
and e* ka = \ {tLka + e Rka) the full Hamiltonian can be 
cast into a compact form 

H = J2 4*4* a k* +J2 V (4«fc + h.c.) + H d , (4) 

ka ka 

with Hd as given in Eq. (1). 

A. Ferromagnetic leads 

For ferromagnetic materials electron-electron interac- 
tion in the leads give rise to magnetic order and spin- 
dependent DOS, p r |(cj) ^ pr^uo). Magnetic order of 
typical band ferromagnets like Fe, Co, and Ni is mainly 
related to electron correlation effects in the relatively nar- 
row 3d sub-bands, which only weakly hybridize with 4s 
and 4p bands 20 . We can assume that due to a strong spa- 
tial confinement of d electron orbitals, the contribution of 
electrons from d sub-bands to transport across the tunnel 
barrier can be neglected 21 . In such a situation the system 
can be modeled by noninteracting 22 s electrons, which 
are spin polarized due to the exchange interaction with 
uncompensated magnetic moments of the completely lo- 
calized d electrons. In mean-field approximation one can 
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model this exchange interaction as an effective molecu- 
lar field, which removes spin degeneracy in the system 
of noninteracting conducting electrons, leading to a spin- 
dependent DOS. 

Parallel and anti-parallel leads ' magnetization. - In ex- 
periments very frequently the electronic transport mea- 
surements are performed for two configurations of the 
leads' magnetization direction 10 : the parallel and an- 
tiparallel alignment. By comparison of electric current 
for these two configurations one can calculate the tunnel- 
ing magnetoresistance (TMR) important parameter for 
the application of the magnetic tunnel junction 12 . 

In this paper we restrict the leads' magnetization di- 
rection to be cither (i) parallel, i.e. the left and right 
lead have the DOS /9l<t(^) and Pr ct (u;) respectively, so 
the total DOS corresponding to the total dispersion e* ka 
(Vfc e [-.Do; -Do]) is given by p a (u>) = Pl<tM + /jr ct (w). 
(ii) antiparallel, the magnetization direction of one of 
the leads (let us consider the right one) is reverted so 
Pr<xM — » Prs(w), where a =| (f) if a =| (J.), so the to- 
tal DOS is described then by p„(co) = Pl„(u)) + p^(u)). 
Here D labels the full (generalized) bandwidth of the 
conduction band (further details can be found below). 
Effects related to leads with non-collinear leads' magne- 
tization are not discussed here 23 . 

For the special case of both leads made of the same 
material in the parallel alignment pha(u) = Pr ct (w) and 
in the antiparallel alignment plct(^) = Prs-(w). There- 
fore for the antiparallel case it gives the total DOS to be 
spin independent p-\(w) = Plt(^) + PRJ.(^) = Plj.(w) + 
/9ri(cj) = P[(u>). In a such situation one can expect 
the usual Kondo effect as for normal (non-ferromagnetic) 
metallic leads, however, the conductance will be dimin- 
ished due to mismatch of the density of states described 
by the prefactor of the integral in Eq. (23) 

B. Different band structures 

In this paper we will consider different type of total 
spin-dependent band structures p a {uj) independently of 
the particular magnetization direction of leads. The par- 
ticular magnetization configuration will affect only the 
linear conductance G a due to the DOS mismatch for both 
leads as discussed in detail later. 



1. Flat band 

The simplest situation, where one can account spin 
asymmetry in is a flat band with the energy-independent 
DOS p a {u) — Pa- Then the spin asymmetry can be pa- 
rameterized just by single parameter the spin imbalance 
in the DOS at the Fermi energy uj = 0. For flat bands the 
knowledge of the spin-polarization P of the leads, defined 
as 

PT (0)- Pi (0) 

pt(0) + M0)' u 



equation is sufficient to fully parameterize p CT (w) (sec 
Fig. 2). This particular DOS shape is a special due to the 
fact, that the particle- hole symmetry is conserved in the 
electrodes leading to particular behaviour for the sym- 
metric Anderson model. This type model of the leads 
will be considered in Sec. V. 

2. Stoner splitting 

One can generalize this model and break the particle- 
hole symmetry by taking the Stoner splitting into ac- 
count. A consequence of the conduction electron ferro- 
magnetism is that the spin-dependent bands arc shifted 
relative to each other (see Fig. 1). For a finite value of 
that shift Ao-, the spin-er band ranges between — D + 
A CT < uj < D + Ao- (with the 'original' bandwidth D). 
In the limit A CT = only energies within the interval 
[— D; D] are considered, a scheme usually used in the 
NRG calculations 7 . This relative shift between the two 
spin-dependent bands leads to so called Stoner splitting 
A, defined as 

A = A L - A T (6) 

at the band edges of the conduction band (see Fig. 10). 
The flat band model with a Stoner splitting and conse- 
quences of the particle-hole symmetry braking is consid- 
ered in Sec. VI. 

3. Arbitrary band structure 

Since the spin-splitting of the dot level is determined 
by the coupling to all occupied and unoccupied (hole) 
electronic states in the leads, the shape of the whole band 
plays an important role. Therefore it is reasonable to 
consider arbitrary DOS shape, which cannot be param- 
eterized by a particular set of parameters as the imbal- 
ance between |- and ^-electrons at the Fermi energy or 
a Stoner splitting (see Fig. 1). Therefore in Sec. VII we 
will consider a model with a more complex band struc- 
ture and in the next Sec. Ill we will develop the NRG 
technique for arbitrary band structure. 

III. THE METHOD: NRG FOR ARBITRARY 
SPIN-DEPENDENT DENSITY OF STATES 

In our analysis we take the ferromagnetic nature of 
the noninteracting leads by means of a spin- and energy- 
dependent DOS P<t(w) into account. A general exam- 
ple is given in Fig. 1. To compute the properties of 
the model described above we have extended numeri- 
cal renormalization group (NRG) technique calculation 
to handle spin-dependent density of states. In order to 
understand to what extent the method applied here is 
different from the standard NRG it is adequate to briefly 
review the general concepts of NRG. 
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FIG. 1: (a) Example of an energy- and spin-dependent leads 
DOS P<j{oj) with an additional spin-dependent shift A a . To 
perform the logarithmic discretization a generalized band- 
width Do is defined. Since we allow for bands with energy 
and spin dependence the discretization is performed for each 
spin-component separately, see panel (b). Since Tied does not 
include spin-flip processes an impurity electron of spin er [cir- 
cles in panel (b)] couples to lead electrons of spin a, a =| (j), 
only. Panel (b) also illustrates that impurity electrons cou- 
ple to leads' electron and hole sates with arbitrary energy ui, 
M < D . 



The NRG technique was invented by Wilson in the 70's 
to solve the Kondo problem 14 - later it was extended to 
handle other quantum impurity models as well 24-26,28 . In 
his original work Wilson considered a spin-independent 
flat density of states for the conduction electrons. Closely 
following Refs[ 17,27] -where the mapping for the case 
of energy-dependent DOS was given- we generalize that 
procedure for the case of the Hamiltonian given in Eq. (4) 
which contains leads with an energy- and spin-dependent 
DOS. 

It is convenient to bring the Hamiltonian given by 
Eq. (4), into a continuous representation 27 before the 
generalized mapping is started. The replacement of the 
discrete fermionic operators by continuous ones, ctka — > 
ow, translates the lead and the tunneling part of the 
Hamiltonian into 

Hi = ^2 / dw g a (u) al CT aw (7) 

Hu = V / diuh a {u) (4aw + h.c.) . (8) 

As shown in Ref. 27 the hereby defined generalized dis- 
persion go-(uj) and hybridization h a {uj) functions have to 
satisfy the relation 



[K {g- x {u))] 2 = p„(u) [V^f , (9) 



where g~ l (u>) is the inverse of g a (uj), what ensures that 
the action on the impurity site is identical both in the 
discrete and the continuous representation. Note that 
there are many possibilities to satisfy Eq. (9). 



The key idea of Wilson's NRG 14 is a logarithmic dis- 
cretization of the conduction band, by introducing a dis- 
cretization parameter A, which defines energy intervals 
] - D A- n ;-D A-' 1 - 1 } and [AjA-"" 1 ; L> A~™[ in the 
conduction band (n € No). Within the n-th interval 



of width dL = A"" ( 1 - A" 



a Fourier expansion of 



the lead operators with fundamental frequency 

Q n = 2tt / d n is defined 



l e ±in„ pw jf A -(«+i) <± UJ< \- 
else 



Here the subscripts n andp (€ Z) label the corresponding 
interval and the harmonic index, respectively, while the 
superscript marks positive (+) or negative (— ) intervals, 
respectively. 

The above defined Fourier series now allows one to re- 
place the continuous fermionic conduction band opera- 
tors a ua by discrete ones a npa {b npa ) of harmonic index 
p and spin a acting on the n-th positive (negative) inter- 
val only 



(10) 



Impurity electrons couple only to the p = mode of 
the lead operators, given the energy-dependent general- 
ized hybridization h a (uj) is replaced by a constant hy- 
bridization, hcr(u>) — > h+ a for u> > (or h~ a for lo < 0, 
respectively). Obviously, the particular choice of con- 
stant hybridization h^ a demands the generalized disper- 
sion g„(u>) to be adjusted accordingly, such that Eq. (9) 
remains valid. Details of this procedure can be found in 
Appendix A. Since we adopt this strategy the harmonic 
index p (the impurity couples only to lead operators of 
harmonic index p = 0) will be dropped below. 
Defining a fermionic operator 



with £ 0(T = J2n 




fine 



(11) 



J j T a (uj)duj 



and 



the coefficients 7^. as given in Appendix A 29 , reveals 
that the impurity effectively couples to a single fermionic 
degree of freedom only, the zeroth site of the Wilson chain 
[for further details see Eq. (A2)]. Therefore the tunneling 
part of H. can be written in a compact form as 



H u = 




h.c. 



(12) 



The final step in the NRG-procedure is the transforma- 
tion of the conduction band Tig into the form of a linear 
chain. This goal is achieved via the tridiagonalization 
procedure developed by Lanczos 30 



Hi 



[Enafnafncr^tncrifnafn+la^fn+icrfncr)]- (13) 
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In general the on-site energies e na and hopping matrix el- 
ements t nt7 along the Wilson-chain need to be determined 
numerically. Besides the matrix elements e na and t n(T 
cocmcients u nma and v nmr7 , which define the fermionic 
operators f na 

oo 

) (14) 

m=0 

already used in Eq. (13), need to be determined. One 
immediately anticipates from Eq. (11) 

U 0m <j = Imcr/V&a J v 0mcr = l ma I ■ (15) 

Equations which determine the matrix elements e na and 
t na and the coefficients u nm a and v nm a are given in Ap- 
pendix B. Note that the on-site energies e niJ vanish in 
the presence of particle-hole symmetry in the leads. 

To summarize: Hamiltonians as the one given in 
Eq. (4) can be cast into the form of a linear chain Hlc 

Hlc =n d + vWtt + foM (i6) 

cr 

oo 

+ y ' [^nerfnerfna + tncr{fncrfn+la + fn+lafncr)] > 
an— 

even though one is dealing with energy and spin- 
dependent leads. In general, however, this involves nu- 
merical determination of the matrix elements e na and 
t na 17 in contrast to Krishna-murthy 24 et al. (who con- 
sidered flat bands) no closed analytical expression for 
those matrix elements is known. 

Eq. (16) nicely illustrates the strength of the NRG- 
procedure: As a consequence of the energy separation 
guaranteed by the logarithmic discretization, the hop- 
ping rate along the chain decreases as t n <~ hr 11 / 2 (the 
on-site energies decays even faster) which allows us to di- 
agonalize the chain Hamiltonian iteratively and in every 
iteration to keep the states with the lowest lying energy 
eigenvalues as the most relevant ones. This very fact un- 
derlines that this method does not rely on any assump- 
tions concerning leading order divergences. 31 



IV. PERTURB ATIVE SCALING ANALYSIS 

We can understand the spin-splitting of the spectral 
function using Haldane's scaling approach 16 , where quan- 
tum charge fluctuations are integrated out. The behavior 
discussed in this paper can be explained as an effect of 
spin-dependent quantum charge fluctuations, which lead 
to a spin-dependent rcnormalization of the dot's level 
position edo- and a spin-dependent level broadening Y a , 
which in turn induce spin splitting of the dot level and the 
Kondo resonance. Within this approach a spin-splitting 
of the local dot level, Aej = Se^ — Se^i + B, which de- 
pends on the full band structure of the leads is obtained 13 



where 

it J I to - 6da eds + U - lo J 

Note that the splitting is not only determined by the 
leads spin polarization P, i. e. the splitting is not only a 
property of the Fermi surface. Eq. (17) is the key equa- 
tion to explain the physics of the (spin-dependent) split- 
ting of the local level e& a . This equation explains the 
spin-dependent occupation and consequently the split- 
ting of the spectral function of a dot that is contacted 
to leads with a particular band-structure. The first term 
in the curly brackets corresponds to electron-like pro- 
cesses, namely charge fluctuations between a single oc- 
cupied state | cr) and the empty |0) one, and the second 
term to hole-like processes, namely charge fluctuations 
between the states | cr) and |2). The amplitude of the 
charge fluctuations is proportional to T, which for T^>T 
determines the width of dot levels observed in transport. 

The exchange field given by Eq. (17) gives rise to pre- 
cession of an accumulated spin on the quantum dot at- 
tached to leads with non-collincar leads' magnetization 
are not discussed here 23 . 



A. Flat band 

Eq. (17) predicts that even for systems with spin- 
asymmetric bands p-\(u) ^ Pi(^), the integral can give 
Ae = 0, which corresponds to a situation where the 
renormalization of e^a due to electron-like processes are 
compensated by hole-like processes. An example is 
a system consisting of particle-hole symmetric bands, 
p a {uj) — p a (—un), where no splitting of the Kondo res- 
onance (Aed = 0) for the symmetric point, ed = —U/2, 
appears. 

For a flat band p a {uj) = p a , Eq. (17) can be in- 
tegrated analytically. For D » U, |ea| one finds: 
Ae ~ (Pr/7r)Re[0(e d ) - <j>{U + e d )], where <j>{x) = 
^(i + ix/2irT) and ^(x) denotes the digamma function. 
For T = 0, the spin-splitting is given by 

showing a logarithmic divergence for £d — > or U + £d — ► 
0. 



B. Stoner splitting 

For real systems p-h symmetric bands cannot be as- 
sumed, however the compensation Ae = using a proper 
tuning of the gate voltage ed is still possible. We can ana- 
lyze also the effect of the Stoner splitting by considering a 
flat-band structure in Fig. 2 using the value of the Stoner 
splitting A = 0.2(-§-)D . Then also from Eq. (17) we 
can expect an additional spin-splitting of the dot level 
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e = o 



B/T=-0.1 



FIG. 2: The leads density of states of a flat [p s (u) = p a ] 
but spin-dependent (p-[ ^ pi) band for spin polarization P = 
0.2. The dark region marks the filled states below the Fermi 
energy. Since Af = Aj, =0 the generalized bandwidth Do = 
D. 



induced by the presence of the Stoner field in the leads 
even for spin polarization P — given by 



Ae 



(St) 



2tt 



■In 



(-e d + Dp - A)(U + e A + D -A) 
(-e d + D )(U + e d + D ) 



G19) 



which for the symmetric level position, ed = —U/2, it 
leads to 



Aef } 



U/2 + e d + D -A 
U/2 + e d + D 



(20) 



Eq. (19) also shows that the characteristic energy scale 
of the spin-splitting is given by T rather than by the 
Stoner splitting A (A 3> T), since the states far from 
the Fermi surface enter Eq. (17) only with a logarithmic 
weight. The difference can be as large as three orders 
of the magnitude, so in metallic ferromagnet the Stoner 
splitting energy is of the order A ~ 1 eV but still the 
effective molecular filed Ae^ St ^ generated by it is a small 
fraction of T - of order of 1 meV so comparable with the 
Kondo energy scale for molecular single-electron transis- 
tors. 10 However, the Stoner splitting introduces a strong 
p-h asymmetry, so it can influence the character of gate 
voltage dependence significantly 

In the next Section we analyze the effect of different 
type of the band structure using numerical renormaliza- 
tion group technique and compare it to that obtained in 
this Section by scaling procedure. 



V. FLAT BANDS WITHOUT STONER 
SPLITTING 

We start our analysis by considering normalized flat 
bands without the Stoner splitting (i. e. Dq = D), as 
sketched in Fig. 2, with finite spin-polarization P ^ 
[defined via Eq. (5)]. 32 In this particular case - of spin- 
(but not energy-) dependent coupling T a = irp a V 2 - the 
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FIG. 3: (a)-(d): Spin-dependent occupation n CT of the local 
level as a function of the leads' spin polarization P for = 
-U/3 [(a) and (c)] and e d = -2(7/3 [(b) and (d)] (related by 
the particle-hole symmetry) for B = (left column) and B = 
—0.1 r (right column). For finite spin-polarization the 
condition n-\ = rij can only be obtained by an appropriately 
tuned magnetic field, as shown in (c) and (d). Parameters: 
U = 0.12 Do, T = (7/6. 



coupling r a can be parametrized via P, IVn = |r(l±P) 
[here + (— ) corresponds to spin j (J.)], where T = F^+T^, 
Leads with a DOS as the one analyzed in this section have 
for instance been studied in Ref. [7,8]. 

Since r CT = irpa-V 2 the spin-dependence of T a can 

be absorbed by replacing V — > V a = V\l h(l ± P) in 



Eq. (4), i. e. a spin-dependent hopping matrix element, 
while treating the leads as unpolarized ones {p a — > p). 
This procedure has the particular advantage that the 
'standard' NRG procedure 24 can be applied meaning that 
the on-site energies (tunneling matrix elements), defined 
in Eq. (13), along the Wilson chain e na vanish, while 
t„ CT 's turn out to be spin-independent. Therefore it does 
not involve the solution of the tedious equations given in 
Appendix B. 



A. Spin and charge state 

Consequently we start our numerical analysis by com- 
puting the spin- resolved dot occupation n a = (n a ), which 
is a static property. Fig. 3 shows the spin-resolved impu- 
rity occupation as a function of the spin-polarization P of 
the leads. Fig. 3(a) and (c) corresponds to a gate- voltage 
of ed = — U/3 (where the total occupation of the system 
n-^+ni < 1) whereas the second line to ed = — 2Z7 /3 (with 
n-f +«x > 1). The total occupation of the system nj 
decreases (increases) for e d — —U/3 (e d = —2U/3) when 
the spin-polarization of the leads is finite P ^ 0. Note 
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FIG. 4: Total spectral functions j4(u) = A ct (oj) for var- 
ious values of the spin polarization P. (a) For td = — U/3 
an increase in P results in a splitting and suppression of the 
Kondo resonance (see inset). The effect of finite P on the 
Hubbard peaks is less significant. The spin-resolved spectral 
function A a (u), shown in (b), reveals that A^(uj) and A^iw) 
differ significantly from each other for P / 0. (c) Spectral 
function for the same values of P as in (a) but for ed = —2/3U 
[due to particle-hole symmetry the results are mirrored as 
compare to (a), however with inverted spins]. Parameters: 
U = 0.12A>, B = 0, T = 17/6. 



that both situations, ed = — U/3 and ed = — 2U/3), are 
symmetric in respect of changing particle into hole states 
and vice versa, which is possible only for the leads with 
particle- hole symmetry. For the gate voltage ed = — U/2 
there is a particle-hole symmetry in the whole system 
leads with a DOS which does share this symmetry even 
in the presence of spin asymmetry. Note, whereas a fi- 
nite spin polarization leads to a decrease in + n± for 
ed > —U/2, cf. Fig. 3(a), it results in an increase in 
n-f + ni for ed < —U/2, cf. Fig. 3(b). Obviously, for 
P = and in absence of an external magnetic field the 
impurity does not have a preferred occupation n-\ = n^. 
Any finite value of P violates this relation: u-f > rij, for 
P > and e d > —U/2 (since <5e dT — T± < T T ~ Se^). 
For ed < —U/2, on the other hand, the opposite behavior 
is found. 

The effect on the impurity of a finite leads polarization, 
namely to prefer a certain spin species, can be compen- 
sated by a locally applied magnetic field, as shown in 
Figs. 3(c) and (d). For e = -U/3 and B/T = -0.1, 
cf. Fig. 3(c), the impurity is not occupied by a preferred 
spin species, n-f = n^, for P ~ 0.5. Due to particle- 
hole symmetry the same magnetic field absorbs a lead 

polarization of P 0.5 for e = -2U/3, cf. Fig. 3(d). 

The magnetic field which restores the condition n-f = , 
henceforth denote as compensation field B comp (P), will 
be of particular interest below. 



• e d =-U/3 

+ £j=-7U/12 




FIG. 5: (a) Compensation field B conl p(P) for different val- 
ues of ed- For flat bands B C om P depends linearly on the leads 
polarization. At the point where there is particle-hole sym- 
metry, so for gate voltage (e d = —U/2), B coulp (P) — for any 
value of P, i. e. the spectral function is not split for any value 
of P even though B — 0. (b) Spectral function for various 
values of P for B = B comp . Note the sharper resonance in 
the spectral function, i. e. a reduced Kondo temperature 7k. 



B. Single-particle spectral function 

Using the NRG technique we can access to the spin- 
resolved single-particle spectral density A a (uj, T, B, P) = 
— —Im.Qd' cr (u)) for arbitrary temperature T, magnetic field 
B, and spin polarization P, where G^(u>) denotes a re- 
tarded Green function. We can relate the asymmetry in 
the occupancy, rtj ^ n^, to the occurrence of charge fluc- 
tuations in the dot and broadening and shifts the posi- 
tion of the energy levels (for both spin up and down). For 
F^0, the charge fluctuations and hence level shifts and 
level occupations become spin- dependent, causing the dot 
level to split 6 and the dot magnetization n-\ — n± to be 
finite. As a result, the Kondo resonance is also spin-split 
and suppressed Fig. 4(b), similarly to the effect of an 
applied magnetic field 33 . This means that Kondo corre- 
lations are reduced or even completely suppressed in the 
presence of ferromagnetic leads. 

Note the asymmetry in the spectral function for P / 
which stems from the spin-dependent hybridization T a . 
In Fig. 4(b), where the spin-resolved spectral function 
is plot, reveals the origin of the asymmetry around the 
Fermi energy of the spectral function. The spectral func- 
tion obtained for polarized leads has to be contrasted to 
that of a dot asserted to a local magnetic field, where a 
nearly symmetric (perfect symmetry is only for the sym- 
metric Anderson model) suppression and splitting of the 
Kondo resonance (around the Fermi energy) appears. 33 
In Fig. 4(c) the gate voltage (ed = —2U/3 ) is chosen 
such that it is particle-hole symmetric to the case shown 
in Fig. 4(a). Due to the opposite particle-hole symme- 
try, the obtained spectral function in Fig. 4(c) is nothing 
else but the spectral function shown in Fig. 4(a) mirrored 
around the Fermi energy. 

Fig. 3(c) and (d) showed that a finite magnetic field B 
can be used to recover the condition n-\ — nj_. Indeed, for 
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that one can easily identify whether the Kondo resonance 
is fully recovered from the positions of A a (u>) relative to 
each other. Since a dip in the total spectral function 
A(u>) = Aa(<jS) is not present for a modest shift of 
A a (tjS) w. r. t. each other, we identify an unsplit Kondo 
resonance with perfectly aligned spin-resolved spectral 
functions henceforth. 

Clearly one can identify from Fig. 6 that the spectral 
function is split for any value ed ^ —U/2. This splitting 
changes its sign at = —U/2. For the particular case 
£d = —U/2, cf. Fig. 6(b), the spectral function reaches 
the unitary limit. For this gate voltage the spin-resolved 
spectral functions A a (uj) have a different height as pre- 
dicted from the Friedel sum rule Eq. (21). 



Friedel sum rule 



FIG. 6: Spin-dependent spectral function A a (u>) for vari- 
ous values of £d, fixed P = 0.2, and B — [(blue) dashed: 
Af(u>), (green) long-dashed: A-\(oj) and (red) dotted A(u) — 
^2 A a (u>)]. The splitting between A-\ and Aj. changes its 
sign at £d = —U/2. The spectral function A(uj) is plot for 
several values of P for ed = —U/3 in Fig. 4 (a). The splitting 
of the spectral function A(uj) depends as well on ed as on P. 
Parameters: U = 0.12D , T = (7/6. 



any lead polarization P a compensation field B comp (P) 
exists at which the impurity is not preferably occupied 
by a particular spin species. To a good approximation 
B comp (P) has to be chosen such that the induced spin- 
splitting of the local level is compensated 37 . One con- 
sequently expects a linear P-dependence of 5 comp from 
Eq. (18). Fig. 5(a) shows the P-dependence of P C omp for 
various values of ed obtained via NRG. The numerically 
found behavior can be explained through Eq. (18). It 
confirms that the slope of B comp is negative (positive) 
for e d > -U/2 (e d < -U/2) and that B comp (P) = for 
£d = -U/2. 

The fact that this occurs simultaneously with the dis- 
appearance of the Kondo resonance splitting suggests 
that the local spin is fully screened at P C omp- The expec- 
tation of an unsplit Kondo resonance in presence of the 
magnetic field B comp is nicely confirmed by our numerics 
in Fig. 5(b) for £d = —U/3. Moreover, the sharpening of 
the Kondo resonance in the spectral function upon an in- 
crease in P (indicating a decrease in the Kondo tempera- 
ture Xk) can be extracted from this plot. The associated 
binding energy of the singlet (the Kondo temperature 
Xk) is consequently reduced. 

Having demonstrated that the Kondo resonance can 
be fully recovered for P C om P even though P / wc 
show that there is also the possibility to recover the un- 
split Kondo resonance via an appropriately tuned gate- 
voltage. In Fig. 6 we plot the spin-resolved spectral func- 
tion for various values of ed for P = 0.2 and B = 0. Note 



A direct consequence of the spin-splitting of the local 
level and the accompanied n-f n\ for and B = 

can be explained by means of the Friedel sum rule, 3 an 
exact T = relation valid also for arbitrary values of P 
and B. This formula relates the height of the spectral 
function (at the Fermi energy) and the phase shift of the 
electrons scattered by the impurity. As already pointed 
out spin-polarized electrodes induce a splitting the local 
impurity level resulting in a spin-dependent occupation 
of the impurity. Thus the knowledge of the local occupa- 
tion can be used to anticipate (via the Friedel sum rule) 
that the local spectral function becomes spin-split and 
suppressed (similar to the presence of a local magnetic 
field) in the presence of spin-polarized leads. 

According to Friedel sum rule, the height of the spec- 
tral function at the Fermi energy A a (0) are fully deter- 
mined by the level occupation n a 

, sin 2 (7r?i (T ) . . 

tt> a = Trn a , A a = (21) 

7Tl a 

where 4> a being the spin dependent phase shift. This 
implies a split and suppressed Kondo resonance for P / 
in absence of a magnetic field. 

An equal spin occupation, n-\ = n^, can be ob- 
tained only for an appropriate external magnetic field 
B = Pcomp- For the latter, in the local moment regime 
(n rs 1) we have n-f = « 0.5, so that tf>^ = (f>^ w 7r/2, 
which implies that the peaks of A^ and A± are aligned. 
Thus, the Friedel sum rule clarifies why the magnetic field 
-Bcomp at which the splitting of the Kondo resonance dis- 
appears, coincides with that for which rif = n^. Eq. (21) 
indicates a new feature namely that the amplitude of the 
Kondo resonance becomes spin dependent (see also the 
discussion in Sec. VE). 

It is important to point out that Eq. (21) is valid only 
for the leads with particle-hole symmetry (see Ref.35). 
For arbitrary DOS shape it is possible to generalize the 
Friedel sum rule. 
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FIG. 7: The isotropic Kondo effect, accompanied by univer- 
sal scaling, can be recovered even for P/0 and B = B comp 
since both the spin-susceptibility Im^s (a) and the spectral 
function (b) collapse onto a universal curve. This statement 
holds for any value of with some corresponding compen- 
sation fields as given in Fig. 5(a). Parameters: U = 0.12Do, 
e d = -17/3, T = U/6. 



D. Spin-spectral function - Kondo temperature 



FIG. 8: (a) Spin polarization dependence of the Kondo tem- 
perature Tk(P) [obtained by applying an appropriate B = 
-Bcomp (-P, £d)]- The functional dependence of Tk(P) found 
via a poor man's scaling analysis (dotted line), see Eq.(6) of 
Ref. 6, is confirmed by the NRG analysis. The plot confirms 
that a dot tuned to the local moment regime, ~^2, a n a — 1, 
shows the universal P-dependence of Tk(P), for arbitrary val- 
ues of e<j. (b) The ed-dependence of Tk for fixed P = 0.2 
at B = -Bcomp. As long as the impurity remains in the lo- 
cal moment regime Haldane's formula 16 for Tk (dotted line) 
properly describes Tk(ed) even though the leads have a finite 
spin polarization. Parameters: U = 0.12 Do, T = U/6. 



We already stated above that the Kondo temperature 
Tk decreases when the leads polarization P is increased. 
Obviously, Tr has to vanish when we are dealing with 
fully spin-polarized leads \P\ = 1. To obtain Tr(-P) we 
calculated the imaginary part of the quantum dot spin 
spectral function 



xU") = r{m)([Sz(t),s z (o)})} 



(22) 



{T denotes the Fourier transform), see also Fig. 7(a), 
in presence of the appropriate B comp (P, e^) and identi- 
fied the maximum in Im^o with Tk(-P). Fig. 8(a) shows 
Tk(-P) normalized to Tr(P = 0) for different values of 
£d- As shown in Ref.[6,7] the functional dependence of 
Tk(-P) can be nicely explained within the framework 
of poor man's scaling. Note that the decrease in Tr 
is rather weak for a modest value of P. In Fig. 7(b) 
we plot, additionally to Xk(-P), the gate- voltage depen- 
dence of Xk(ed) normalized to 7k(ed = —U/2) (where 
the Kondo temperature is minimal) for fixed P = 0.2. 
In presence of -B comp (P, ed) this functional dependence 
of the Kondo temperature can analytically be described 
via Haldane's formula 16 for the Kondo temperature of 
a single level dot coupled to unpolarized leads, Tr = 
^VUTe^ ed ^ d+u)/ru . This fact is another manifestation 
that indeed the usual Kondo effect can be recovered in 
presence of spin-polarized leads, given the appropriate 
-Bcomp (-P, £d) is applied. 

Having shown that the coexistence of spin-polarization 
in the leads and the Kondo effect is possible (given 
a magnetic field B is tuned appropriately at a given 
gate- voltage £d), we plot the properly rescaled spin- 
susceptibility Imx|, see Fig. 7(a), and spectral function, 
see Fig. 7(b), confirming that the isotropic Kondo ef- 
fect is recovered. The expected collapse of Imxf and 
Y] T a A a {ui) onto a universal curve for temperatures be- 



low T K [note that T K depends on P, see e. g. Fig. 5(b)] 
is confirmed by our numerical results. For energies much 
bigger that Tr, to 3> Tk, universal scaling is lost and the 
curves start to deviate from each other. 



E. Conductance 

The knowledge of the spectral function enables us 
to compute quantities which are experimentally acces- 
sible, such as the linear conductance. To investigate the 
question whether spin-polarized ferromagnetic leads in- 
troduce a spin-dependent current we compute the spin- 
resolved conductance G a , 



Go 



s 2 2r T ,„r 



Lcr-L Rcr 



fl (Tlo- + Tro-) J_ 



duo A G {u)) 



duj 



(23) 



with f(u>) denoting the Fermi function (remember that 
we choose Tlo- = Tro-). The total conductance G is noth- 
ing else but the sum of the spin-resolved conductances 
G = J2 a Go- Based on the Friedel sum rule, Eq. (21), we 
expect the height of the spin-resolved spectral function 
at the Fermi energy ^(0) ~ l/T a ~ 1/(1 ± P). This 
expectation is nicely confirmed by our numerical results 
in Fig. 9(a). As the T = conductance is given by the 
product of A G (0) and T a the spin-resolved conductance 
G(j ~ r <T ^4 (T (0) becomes P independent. The results of 
the NRG calculation are shown in Fig. 9(b) confirming 
this expectation. We conclude that, even though one is 
dealing with spin-polarized leads, it is not possible to 
create a spin current via spin polarized leads the current 
remains spin independent. For the antiparallel alignment 
due to the DOS mismatch, the conductance will be re- 
duced below Go = 2e 2 /h limit. 
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FIG. 9: (a) Height of the spin-resolved spectral function 
A a (to = 0) as a function of P. The dashed line shows the 
1/(1 ± P) dependence as expected from the Friedel sum rule, 
Eq. (21). (b) Spin-resolved conductance G a as obtained from 
Eq. (23). This plot confirms the expectation based on the 
Friedel sum rule that the spin-resolved conductance should 
be independent of the spin species, consequently it serves as 
a consistency check of our numerics. 



VI. FLAT BANDS WITH STONER SPLITTING 

In ferromagnetic leads, additional to a finite lead spin 
polarization P at the Fermi surface, a splitting between 
the spin-| and spin-| bands A CT (Stoner splitting), which 
is effect of the effective exchange magnetic field in the 
leads, might appear. Taking this effect into account we 
do not expect the splitting of the Kondo resonance to dis- 
appear at ed = —U/2 as in Fig. 6(b), rather at a different 
value of ed determined by the shifts A CT . This is due to 
the fact that the Stoner splitting breaks the particle-hole 
symmetry in the electrodes. Consequently we shall con- 
sider flat, spin-polarized leads (of polarization P), which 
are shifted relative to each other by an amount A a in this 
section. Since particle-hole symmetry is lost for this band 
structure in the leads we expect neither the splitting of 
the Kondo resonance, see e. g. Fig. 6, nor the compensa- 
tion field -Bcomp to be symmetric around e c \ = —U/2, in 
contrast to Section V. 

To quantify the statement that a shift A CT introduces 
an finite exchange field and splitting, in this section we 
analyze the particular case Aj = — A| = A, sketched 
in Fig. 10. Note that the leads DOS structure discussed 
in this section 38 already requires the use of the extended 
NRG scheme, explained in Section III. In particular, 
spin-dependent on-site energies along the Wilson chain 
e na , as explained in Section III and Appendix A and 
B, need to be determined to solve for the leads band 
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FIG. 10: The DOS in case of flat bands with the Stoner split- 
ting A and a finite spin polarization P. In general the shifts 
of the | and \ bands are unrelated to each other, A| 7^ — Af. 
In this section, however, we assume A| = — Af. 



0.75 



0.25 



0.75 



0.25 





CQ 




• e d=- 


-U/3 


♦ e d=- 


-U/2 


+ e d=- 


-7U/12 



-0.5 0.5 
P 



FIG. 11: Spin polarization P dependence of the local occupa- 
tion n a for 73 = (a) and B = 0.1 T (b) with e d = -17/3. The 
Stoner-splitting A even in the absence of spin polarization of 
the leads, P = 0, introduces an effective exchange field which 
splits the local level nj > n-f in absence of an external mag- 
netic field B = 0. The spin polarization dependence of 73 comp 
for various values of ed is shown in (c). To a good approx- 
imation the Stoner-splitting introduces a constant exchange 
field (here: -B com p/r ~ 0.061), therefore it does nothing else 
but shifting the B comp -dependence of section V. Parameters: 
U = 0.12 (-^) D (h T = U/6, A x = — A T = 0.10 (£) D . 



structure shown in Fig. 10. 

The studies summarized in Section V revealed that the 
spin-resolved impurity occupation n a is the key quantity 
in the context of a dot contacted to spin-polarized leads. 
Following this finding we plot n a vs. P for leads with 

A = 0.1 D and e d = -C7/3 in Fig. 11. In contrast 
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FIG. 12: Spin-resolved equilibrium spectral function 
A a (u>, T,V = 0) with parameters extracted from Ref. 11. The 
(blue) dashed lines correspond to At(u, T, V = 0), the (green) 
long dashed lines to Aj(w,T, V = 0) and the (red) solid ones 
to the sum of both, i. e. A(u>,T,V — 0). (a)-(e) A gradual 
increase in the temperature for those temperatures used in 
Ref. 11. For comparison we plot the T = spectral func- 
tion A(u), T = 0, V — 0) in panel (a) as well. Note that the 
splitting of the Kondo resonance disappears upon increasing 
T. 



to the previous section the impurity is preferably occu- 
pied by .[-electrons in the absence of a magnetic field, see 
Fig. 11(a) even for P = 0. A finite magnetic field, e. g. 
B/T = 0,1 [Fig. 11(b)], has the same consequence as de- 
scribed before, namely the shift the local impurity levels 
£dcr- The P-dependence of B comp for the particular band 
structure sketched in Fig. 10 is shown in Fig. 11(c); it is 
roughly given by the compensation field for flat bands, 
Fig. 5(a) where the same gate- voltages were used, shifted 
by an effective exchange field Ae d St ' ) generated by the 
band splittings for P = 0. The value of -Bstoncr can be 
obtained via integrating out those band states of energy 
w which lie in the interval D < \oj\ < Dq (see Eq. (20). 
The value of this effective exchange field, which is loga- 
rithmically suppressed (as can be shown by perturbative 
scaling), given by Eq. (20). Inserting the numbers used 

(St) 

in this section we obtain a Ae d /V as 0.060. The numcr- 

(St) 

ical calculation reveals Ae d /V w 0.061 [see Fig. 11(c)], 
i. e. the numerical result agrees reasonably well with the 
result based on scaling analysis (See Sec. IV). 



A. Finite temperature: Comparison with 
experimental data 

In a recent experiment of Nygard 11 et al. an anoma- 
lous splitting of the zero bias anomaly in the conduc- 
tance induced by the Kondo resonance in absence of a 



magnetic field was observed. In this experiment quan- 
tum dot based on single wall carbon nanotubes (SWNT) 
contacted to non-magnetic, spin unpolarized Cr/Au elec- 
trodes were used. Those authors related the observed 
splitting of the Kondo resonance at zero magnetic field 
with the presence of a ferric iron nitrate nanoparticle 39 , 
which due electronic tunnel coupling to the dot intro- 
duces a spin-dependent hybridization. 

We model the setup of Ref. 11 by means of a single- 
level dot tuned to the local moment regime ed CT = —U/2 
[dashed line in Fig. 1(b) of Ref. 11]. Moreover we used 
Tk = 500 mK 40 , which is in roughly the value of Xk ob- 
served in Ref. 11, as an input parameter and extracted 
from this value U/T — 4.3. Since a splitting of the 
Kondo resonance was observed for e^a = —U/2, see e.g. 
Fig. 1(b) of Ref. [11], a particle-hole asymmetry must ex- 
ist to achieve a splitting. We try to model this effect by 
using flat band structure and by introducing the Stoner 
splitting. We find the best agreement between the exper- 
imentally observed splitting of the Kondo resonance for 
A/U = 2.5. The presence of the nanoparticle introduces 
a spin-dependent hybridization which we parameterize 
via the leads spin polarization P. 

The comparison of the dl/dV vs. V characteristics 
of Ref. [11] with the theoretical result turns out to be 
rather complicated as it involves the nonequilibrium spec- 
tral function A(u>, T, V) which itself is not known how to 
compute accurately 41,42 It is possible to compare quali- 
tatively the splitting in the non-equilibrium conductance 
obtained in the experiment to the single particle spec- 
tral function. Due to the splitting in the particle spectral 
function for the dot, that is strongly related with the low- 
bias conductance measurements we find close similarity 
between both of them. When we approximate the exper- 
imental results of dl/dV with A(ui, T, V = 0) we achieve 
the best agreement between theory and experiment for 
P = 0.1. For these parameters (e d(7 = -U/2, U/T = 4.3, 
A/U = 2.5 and P = 0.1) we compute the temperature 
dependent equilibrium spectral function A(u, T,V = 0) 
presented in Fig. 12. 

The numerical findings shown in Fig. 12 qualitatively 
confirm the behavior found in the experiment of Nygard 
and collaborators, namely vanishing splitting in the 
dl / dV curves upon with increasing temperature. In the 
experiment, however, the differential conductance dl / dV 
turns out to be asymmetric around the Fermi energy [see 
Fig. 2(b)-(e) in Ref. 11]. We attribute the asymmetry in 
the dl / dV curve with by the asymmetry in the coupling 
between left and right leads. 

Finally we want to compare the temperature depen- 
dence of the linear conductance G found in Ref. [11] (see 
Fig. 2(a) of Ref. [11]), with the one based on the model 
described above. In absence of a magnetic field B = 
and source-drain voltage V = a plateau in G was 
found around T w 500 mK in the experiment. Note 
that the calculation of G does not involve the nonequi- 
librium spectral function so we can calculate it exactly. 
In Fig. 13(a) we show the theoretical curve of the spin- 
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FIG. 13: Linear conductance as a function of temperature (a) 
for the model studied in this section and (b) for the case of 
a dot coupled to normal spin unpolarized leads (P = and 
A = 0) but in presence of a magnetic field B — 0.0217. In 
both cases a plateau in the linear conductance around 7k can 
be observed. Whereas the spin-resolved linear conductance 
is degenerate in case (b) (since it is the symmetric Anderson 
model with particle-hole symmetry), this is not the case for 
(a). Parameters as in Fig. 12, besides P — A = in case (b). 




-1 



-1 



FIG. 14: The parabolic density of states (typical for s electron 
band) given by Eq. (24) with the same Stoner splitting A = 
0.3D but with an additional spin asymmetry Q defined in the 
text. Here Q = 0.0 (a), 0.1 (b), and 0.3 (c). 



For comparison we also plot the linear conductance of 
a dot contacted to normal spin unpolarized leads (i. e. flat 
bands with P = which are not shifted relative to each 
other) for increasing temperature, however in presence of 
a finite magnetic field B — 0.02 U in Fig. 13(b). Such a 
magnetic field might e. g. be inserted in the system via 
the presence of the ferromagnetic nanoparticles. Also 
such a scenario leads to a plateau in G. The value of 
the associated magnetic field, however, is rather big B = 
0.02 U corresponding to B s» 1 T [here we assumed g = 
2.0 for the conduction electrons in the nanotube (see also 
Ref. 11)]. According to Ref. 11 the latter scenario can 
therefore not explain the observed behavior, since one 
can argue that such a large magnetic field is hardly to 
believe can exist in this experiment (see in particular the 
footnote 13 of Ref. [11]). 

Consequently, the presence of ferromagnetic leads 
might explain the experimental observation. We can, 
however, not exclude other possible mechanisms - such 
as the existence of an effective magnetic field - which 
might lead to the observed behavior, however, only form 
special geometry, where the strong stray magnetic field 
is possible. 



VII. BAND STRUCTURE WITH ARBITRARY 
SHAPE 

Flat bands, studied in the previous sections, are only 
a poor approximation of a realistic band-structure in the 
leads. 43 We complete our study by analyzing the effect of 
an arbitrary leads DOS. Consequently we are considering 
leads with a DOS which is energy- and spin-dependent, 
and additionally contains a Stoner splitting in this sec- 
tion. As outlined in the previous section the 'generalized' 
NRG-formalism introduced in Section III has to be used 
in this case as well. In particular, we study the effect of 
gate voltage variation on the spin-splitting of the local 
level of a quantum dot attached to ferromagnetic leads. 
We show how the gate voltage can control the magnetic 
properties of the dot. A similar proposal, namely to con- 
trol the magnetic interactions via an electric field, was 
recently made in gated structures. 44 

In an analysis of several types of band structures 
we found three typical gate voltage dependences of the 
Kondo resonance. To be more specific, to illustrate this 
behavior we used a square-root shape DOS or parabolic 
band (as for free s-type electrons) 45 with the Stoner split- 
ting A 46 , and some additional spin asymmetry Q 



13a/2 



£>" 3/2 (l + o-Q) yjuj + D + aA/2 (24) 



resolved conductance G a obtained via Eq. (23) for the 
model discussed in this section. In agreement with the 
experiment a plateau in the linear conductance G is found 
around Tk- 



depicted in Fig. 14 [a = 1(— 1) for f {{)]. This example 
turns out to encompass three typical classes mentioned 
above. Note that we restricted the DOS in Eq. (24) u> £ 
[-D — aA/2,D — crA/2] to ensure that we are dealing 
with a normalized DOS. 
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FIG. 15: Color scale plot of the (rescaled) spectral function 
y)„ T a A a (uj): the dashed lines mark the Kondo resonance for 
the leads band-structure as shown in Fig. 14 in absence of an 
external magnetic field. Whereas panel (a)-(c) correspond to 
a quantum dot in absence of an external magnetic field, in 
panel (d)-(f) an external field is applied. 



FIG. 16: Color scale plot of the (rescaled) spectral function 
^ r CT y4 CT (aj): the dashed lines mark the Kondo resonance for 
the leads band-structure as shown in Fig. 14. In contrast to 
Fig. 15 an external field is applied. As explained in Section V 
a local magnetic field shifts the spin-resolved local levels rel- 
ative to each other leading to a change in the gate voltage 
where the compensation (rif = nj.) takes place. 



The three typical gate voltage dependences of the 
Kondo resonance we found: (i) a scenario where the 
splitting of the local level was roughly independent of 
gate voltage [Fig. 15(a)], (ii) a scenario with a strong 
gate voltage dependence of the splitting, however with- 
out compensation [Fig. 15(b)] and (hi) a case with a 
strong gate voltage dependence of the splitting includ- 
ing a particular gate voltage where the splitting vanishes 
[Fig. 15(c)]. Fig. 16(a), (b), and (c) shows the effect of 
an external magnetic field for Fig. 15(a), (b), and (c) re- 
spectively. As already outlined in the previous sections 
an external magnetic field can be used to compensate the 
lead induced spin-splitting (which itself depends on gate 
voltage) and thereby change that gate voltage where the 
full Kondo resonance exists, see Fig. 15(c) and Fig. 16(c). 



The spitting of the spectral function (see Figs. 15 
and 16) agrees very well with the splitting based on Hal- 
dane's scaling approach 16 given by the formula Eq. (17). 
The white dashed line in Figs. 15 and 16 shows prediction 
of Eq. (17), which fit very well to the numerical results 
presented. 

In Fig. 17(a)-(c) we plot the linear conductance 47 G for 
the three band-structures sketched in Fig. 14 as a func- 
tion of gate voltage £d and external magnetic field B. 
The two horizontal maxima in the conductance around 
£d ~ and £d ~ — U correspond to standard Coulomb 
resonances. The increased conductance in between, how- 
ever, is due to the Kondo effect. Note that, in contrast 
to the Kondo effect observed in case of a dot attached to 
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FIG. 17: Color scale plot of the linear conductance G as 
a function of gate voltage e c \ and magnetic field B for the 
band-structures corresponding to Fig. 14(a)-(c) [i. e. (a) cor- 
responds to Q = 0, (b) to Q = 0.1 and (c) to Q = 0.3]. The 
black dashed line indicates cross-sections for which the am- 
plitude is plotted in Fig. 19(b). Note the charging resonances 
around ea w and ed w — U and the tilted resonance (due to 
the Kondo effect) for -U < e A < 0. 



flat and unpolarized leads (where the Kondo resonance 
is a vertical line), the Kondo resonance has a finite slope 
here. This is due to the gate voltage dependence (and 
correspondingly to the compensation field dependence) 
of the the spin splitting of the local level [Eq. (17)]. 

It is interesting to consider how the Kondo resonance 
merge with Coulomb resonances. As one can learn form 
Eq. (18 ) that the splitting Sea for the flat band structure 
without Stoner splitting shows a logarithmic divergence 
for eo — * or U + eo — ► 0. Since any sufficiently smooth 
DOS can be linearized around the Fermi surface, this 
logarithmic divergence occurs quite universally, as can be 
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FIG. 18: Log-linear version of the Fig. 17(c) - color scale plot 
of the linear conductance G as a function of gate voltage e<j 
and magnetic held B for the band-structures corresponding 
to Fig. 14(c). 
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FIG. 19: The spin-resolved impurity occupation n„ for the 
DOS reflecting Fig. 14(a) and Fig. 15(a) is plot in panel (a). 
It nicely illustrates that £d can be employed to tune the av- 
erage impurity magnetization. Since the Kondo resonance is 
only fully recovered for rif = ny the linear conductance G 
shows three peaks in panel (b) for three different situations 
corresponding to black dashed lines from Fig. 17. 



observed in log- linear versions of Fig. 17(c) demonstrated 
in Fig. 18. For finite temperature (T > 0) the logarithmic 
divergence for eo — > or e — > — U is cut off, Ae ~ 
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-(l/7r)Pr[tf(l/2)+ln(27rT/t/)], which is also important 
for temperatures T<Tk. 

The finite slope of the Kondo resonance in the linear 
conductance plot motivates us to suggest an interesting 
possibility here: a dot attached to leads with a DOS as 
described in Eq. (24) and say Q = 0.3, see Fig. 17(c), the 
gate voltage can be used to tune the magnetization of the 
dot. In absence of an external magnetic field B = the 
dot is not occupied by a preferred spin species for a par- 
ticular gate voltage (which is w —U/2 in this case). 
Consequently the dot is preferably occupied with elec- 
trons of spin a for ed > — U/2 and electrons of opposite 
spin a for ed < —U /2. In other words, the average spin on 
a dot contacted to leads with such a DOS, n a , and conse- 
quently its magnetization m = n^ — ni arc tunable via an 
electric field, which in turn affects ed- This observation 
suggests the interesting possibility (similar to Ref. 44) 
to employ an electric field to tune the magnetic proper- 
ties of a dot. To compute the gate voltage dependence 
of the magnetic field, of course, a detailed knowledge of 
the band-structure in the leads is necessary. Fig. 19(a) 
summarizes this discussion. 

Another consequence of the leads band-structure is 
that the linear conductance G does not show a plateau 
when plotted relative to ed- Instead, we observe a three 
peak structure in G [see Fig. 19(b)] due to the partic- 
ular gate voltage at which the full Kondo resonance is 
recovered. 

The theoretical predictions made in this section were 
done for particular leads DOS of Eq. (24). We used this 
particular leads DOS to illustrate the effects on the spin- 
splitting of the local level triggered by the band-structure 
of the leads. To quantitatively compare our predictions 
with experiments, however, a detailed knowledge of the 
leads DOS is necessary which is in general very difficult 
to obtain. 



VIII. CONCLUSION 

We analyzed the effect of the band structure of the fer- 
romagnetic leads on the spin splitting of the local level 
edo- of a quantum dot attached to them. We analyzed 
the effect of a finite spin polarization P in case of flat 
(energy-independent) bands. We found that the a finite 
lead polarization results in a spin splitting, of the local 
level and the Kondo resonance in a single particle spec- 



tral function, which can be compensated by an appropri- 
ately tuned magnetic field. Given the condition nj = 
is recovered the isotropic Kondo effect is recovered even 
though the dot is contacted to leads with a finite spin po- 
larization. Additional to this we identified that a finite 
Stoner splitting introduces an effective exchange field in 
the dot. Finally we explained the consequences of an 
energy- and spin-dependent band-structure on the local 
level. We confirmed that also in this general case the 
Kondo resonance can be recovered. From a methodolog- 
ical point of view we extended the 'standard' NRG pro- 
cedure to treat leads with an energy- and spin-dependent 
DOS. 
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APPENDIX A: TRANSFORMATION OF 
TUNNELING HAMILTONIAN - H ld 

The continuous representation of the tunneling Hamil- 
tonian Hid, Eq. (3), can be rewritten in terms of discrete 
conduction band operators for each spin component sep- 
arately by replacing the continuous conduction band op- 
erators ct ua by discrete ones, as suggested in Eq. (10) 



A-(«+i) 



du h a {u)^+ p {u) 



/ 



duj n _(w) +h.c. 



.(Al) 



As outlined in Section III we replace the energy depen- 
dent generalized hybridization of the n-th logarithmic in- 
terval h a (ui), A _( ™ +1) ± lu < A~™, by a constant 
defined as 



1 else 



Ka = ( £ r A A ? +1> [V^)f if -A"" < w < -A-(«+ ] 

1 else 



Consequently possible integrals in Eq. (Al), such as for p = (due to the Riemann-Lebesgue Lemma) a fi- 
/a-(„ + d h a (w)*+ p (u>) oc Xf-c+D d^npM, g ive onl y 
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nitc contribution. For the particular choice of constant 
hybridization the impurity couples to s-waves (p = 
modes) only; consequently we skip the harmonic index 
p below. With above definitions Eq. (Al) simplifies to 
the following compact form 



Hid 



where 



-y 



4 E ( a na"fna + b nal na ) + h - C - 



(A2) 



In 



A -(n + l) 



(A3) 



Remember that Eq. (9) forces the generalized disper- 
sion 5,7 (e) to be adjusted accordingly for this particular 
choice of h a (u>). It is explained in details in Appendix B. 



APPENDIX B: MAPPING OF CONDUCTION 
BAND ONTO SEMI-INFINITE CHAIN 

Here we outline the important steps to bring the leads 
Hamiltonian He, Eq. (8), into the tridiagonal from in- 
troduced in Eq. (13). The replacement of the continu- 
ous conduction band operators aw by discrete opera- 
tors, cf. Eq. (10), simplifies He significantly. 48 As shown 
in Ref. 27 the particular choice of the generalized hy- 
bridization function h a (uj) — ► h„ a (given in Appendix A) 
results in 

^ = E [^LaW + Qa^naKcr] , (Bl) 



where 



J A - ( n+i) * P<r(e) 
J_A-n & £ Me) 

- r _A-(«+i) , , , 



(B2) 



(B3) 



To achieve the goal to bring He into tridiagonal form the 
tridiagonalization procedure developed by Lanczos 30 is 
used. In general diagonal and off-diagonal matrix ele- 
ments e na and t ncr 49 need to be computed in the course 
of the procedure. These matrix elements can be obtained 
by demanding the following relation (see e.g. Ref. 17) 

E (Gk 4<rGW + Cn* bi„b na ) = (B4) 

na 

= ^ ] ^nafnafna + ^ncr (^.fnafn+la + /rt+lcr/nf^ • 

The spin-dependent coefficients u nmi7 and i> nm(T of the 
single-particle operator [Eq. (15)] (that acts on the 



n-th site of the Wilson chain) f n<7 , given by the 
Ansatz Eq. (14), need to be determined recursively for 
each spin component separately. Inverting the Ansatz 
Eq. (14) enables one to express the discrete conduction 



m — u mna J ma 



f ma and b n 



band operators as a na = 

Y^m=o v mno- fm.cn respectively. When we insert a na and 
b na in the l.h.s. of Eq. (B4) and compare corresponding 
j na operators on both sides of this equation we obtain 

oo 

{Cms u nma0^ m a + Cms v nmab ma ) = 



m=0 



— e na.fna + ^nafn+la + ^n-la f\-l a ■ (B5) 

In particular, Eq. (B5) for n = gives the relation for 
the operators foa yields 



E 

m=0 



C + 7+ + 
rr — a n 



+ 



(ma'Ymo 



£,0a 



bt 



(B6) 



where the corresponding values of uo ma - and vo ma 
[Eq. (15)] have already been inserted. Since the oper- 
ators f na obey Fermi-statistics, {/ TO ,/* V } = <W<W, 
the anticommutator of the r.h.s. of Eq. (B6) with /o CT 

yields {e afoa + t afl a , ha} = e 0a . The corresponding 
anticommutator of the l.h.s. of Eq. (B6) with foa finally 
results in 

£ ° CT = T^E fc™ fc) 2 + Cma (7ma) 2 l ■ (B7) 

The initial hopping matrix element to a is readily obtained 
from the anticommutator {eocr/cL + t afl a , eoaha + 
W/icr} = (eotr) 2 + (W) 2 - As e (T is known, t 0CT can be 
easily obtained from Eq. (B6) 



(t0af = £ [(CL) 2 (7+ CT ) 2 + (Cmaf (imaY 



E [C- + C mff (7 mCT ) 2 ] / MB8) 

The knowledge of e 0(T , £o<x, and / cr now enables us to 
extract the coefficients of fi a from Eq. (B6) 



7^ 



'^Oatoa 
~1ma 
'(.Oatda 



((ma - e 0a) , 
(Cma - e 0a) ■ 



(B9) 



A generalization of the argumentation above finally al- 
lows one to obtain the spin-dependent on-site energies 
e na and hopping matrix elements t na for the n-th site of 
the Wilson chain 

e n<r = ^ ] [(fnmu) Cma { v nma) Cma\ i (BIO) 
m 

(incr) = ^(w nmcr ) (Cmcr) + Kmi) (Cncr) 

m 

- (t„-l CT ) 2 - (enaf , (Bll) 
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where the involved coefficients of the single-particle op- 
erator f n +i a , are given as 

u n-\-lmcr — T [(Cmcr "~ e ncr) tinnier ^n — la^n— lma\ ; 
^n+lmcr — 7 [(Cmcr — e ncr) ^nma ~ ^n— lcr^n— lmcr] • 

(B12) 

Note that it is crucial to determine the coefficients 



Unma and v nma even though, only the matrix elements 
along the Wilson chain (e ni7 and t niJ ) arc finally required. 
The numerical solution of the above mentioned equations 
is rather tricky: since the band energies are exponentially 
decaying in the course of the iteration one has to use re- 
liable numerical routines (i.e. arbitrary-precision Fortran 
routines 27 ) to solve for the coefficients u nma and v nma 
and the matrix elements t nrT or e na , respectively. 
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